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Abstract 

Here we show how to generate a dark two-mode squeezed state of a trapped ion, employing 
a three-level ion in a V configuration with a strong decay of the excited states. The degree of 
squeezing can be manipulated by choosing the intensity of the driving fields. Our scheme is robust 
against the usual dissipation mechanism and could be implemented with present-day technology. 
The validity of the approximations employed in this work was tested by numerical calculations, 
which agreed completely with the analytical solutions. 
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The recent experimental advances in Quantum Optics, specially in the domain of trapped 
ions, have allowed fundamental features of quantum mechanics, such as geometric phases 
[l| and Bell inequalities 2 L to be investigated, as well as offering potential applications in 
quantum computation 3, 4] and teleportation processes [5!]. With the advent of quantum 
information theory, the generation of entangled states has became essential for the imple- 
mentation of quantum communication protocols [3J] and to improve our understanding of 
this non-local character of the quantum theory 6]. In particular, the two-mode vacuum 
squeezed (TMVS) state, i.e., the original Einstein- Podolsky- Rosen (EPR) state [7], has at- 
tracted much attention because it can show a high degree of entanglement and can be 
useful for teleportation of continuous variable states |9|]. Success in generating the TMVS 
state has been reported in the running wave domain, with a parametric down conversion 
process However, the experimental generation of this state in the cavity quantum elec- 
trodynamics (QED) or trapped ion domains has not been achieved so far, mainly due to 
the sensitivity of quantum states to system-environment interaction. In the cavity QED 



context, several theoretical schemes with three-level atoms lid . I 111 . [12fl , or even two-level 
atoms, where the sideband transition is used [13|], have been elaborated for the generation 
of the TMVS state. Also in the trapped ion domain we find some schemes which allow 
the generation of this entangled state through the manipulation of laser fields 14J. How- 
ever, none of the schemes cited above take into account the system-environment interaction, 
which degrades the quantum states so that, in general, the fidelity of the generated states 
decays quickly. In this scenario, reservoir engineering appears to offer a possible way round 
this problem and can generate robust non-classical states of the radiation field or of the 
ionic motion. For example, using the atomic decay of the internal levels of a single ion, in 



Ref. [15[ the authors showed how to construct a reservoir able to lead the motion of the 
ion to an squeezed state asymptotically. Similar schemes have been employed to protect 
various superpositions of coherent states Q, and in ref. [ljl the authors showed 

how to protect any one-dimensional motional state of an ion against decoherence. Also in 
this context of reservoir engineering, in Ref. 2o| we find an effective master equation that, 
in the stationary state, filters specific number states of the vibrational motion of a trapped 
ion. On the other hand, reservoir engineering for multi-mode states has been addressed 
only recently, and there are few theoretical schemes so far. We can cite, for example, in the 
trapped ion domain, theoretical schemes for the preparation of a pair coherent state 21[ and 
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pair cat states 



22|, SU(1,1) inte 



ligent states [23J and dark SU(2) states of a trapped ion 
241 ] . Recently, Parkins et al. 25] has proposed a scheme for the unconditional generation 
of a two-mode squeezed state of two separated atomic ensembles. A similar scheme was 
emp loyed for the generation of the TMVS state for the motion of two ions in different traps 
261 ] or even a single ion in a two-dimension trap 27|] inside an optical cavity. In Ref. 28j] it 
was shown theoretically how to generate this entangled state using an atomic reservoir for 



a two-mode cavity. In Ref. [29|, the authors showed how a beam splitter operation may be 
produced in a single ion in two-dimension trap. After generating a robust squeezed state of 
a single mode of a trapped ion |15|], this effective interaction could be directly employed to 
generate a TMVS state, but in this case the TMVS state would not be the steady state of 
the system. 

In this communication we report a simple feasible scheme for the unconditional generation 
of the TMVS state in a single trapped ion. For this purpose we have employed a two- 
dimensional harmonic motion (on the x and y axes) of the center of mass of a single ion in 
a V configuration (see Fig. 1). The excited states, |1) and |2), are coupled to the ground 
state |0) through classical fields (propagating along the x and y axes). When the decay of 
the excited electronic states is stronger than the effective coupling between the vibrational 
and the internal ionic states, the steady state of the vibrational modes, for convenient 
choices of the intensity and frequency of the classical fields, turns out to be exactly the 
TMVS state. Even in the presence of a thermal reservoir for the ionic motion, we show 
that the generated TMVS state is almost exactly the desired one. Such a scheme, based 
on reservoir engineering, is robust against dissipative effects of the vibrational modes and 
could be used to investigate experimentally the entanglement properties of this state. The 
basic level configuration needed for the implementation of our scheme is sketched in Fig. 1. 
We consider an ion with mass m in a two-dimensional trap, driven by four classical fields, 
two of them along the x axis and the other two along the y axis, with complex amplitudes 
Qja = \&ja \ & ltpia (\tya\ being the Rabi frequency and <pj a the phase of the classical fields), 
frequencies Uj Q , and wave numbers ki a , j = 1,2 and a = x,y. The total Hamiltonian of the 
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system is given by E = E + V(t), with 

Hq = fajj\0\\ + hw 2 a 22 + hv x a)a + hiSytfb, 

v(t) = h [n lx e ikixX ~ iuJixt + n ly e ikl * v - iul » t ] <7 10 

+ h [Q 2x e ik2xX ~ iuJ2xt + Qzye^yv-^] a 20 + h.c, 



(la) 



(lb) 



where oj\ and u 2 stand for the atomic transition frequencies between the states |0) < — > |1) 
and |0) < — > |2) respectively, <Ji m = \l) (m\, l,m — 0,1,2, are the atomic operators, a (6) 
and a) (&') are the annihilation and creation operators of the vibrational mode in the x (y) 
axis, with frequency v x {v y ), and h.c stands for Hermitian conjugate. In the Lamb-Dicke 
limit, i. e., r)j a <C 1, where rjj a = kj a J 2m ^ , j — 1, 2 and a = x, y, the above Hamiltonian 



can be written in the interaction picture as 



E 1 = htt lx e- iSl * t [1 + irj lx (ae- iVxt + oV"**)] <r w 
+ m ly e- iSlyt [1 + irii y (be'™** + feV^*)] a 10 
+ nn 2x e- i5 ^ [1 + irj 2x (ae-*"* + aV 1 "*')] a 20 
+ n 2y e- i6 ^ [I + i V2y (be-**** + tie**')] a 20 + /i.e., 



(2) 



with <5 iQ = - LOi. Supposing 5 lx = -8 2x = -v x , 5 ly = -5 2y = v y , \5 ia \ > \r) ia Q ia \, and 
applying a rotating-wave approximation, the effective Hamiltonian becomes 



E! = h\ lx <a + 



A 



ly 



6 f | cr 10 + h\ 2y | b 



+ ^a 1 " ^ a 20 + /i.e. 



(3) 
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where we have defined Xj a = ir)j a Qi a . As in Ref. 25], we can apply a unitary transformation 
P = Sab(0pSab(£), with S a fc(£) = exp (£*afe — ^aV) and £ = e J<? V, this last being the two- 
mode squeezing operator (r stands for the squeezing factor and the angle of squeezing), 
to obtain the transformed Hamiltonian 



Ei = Sl^EjSab^) = h\ a aa w + h~\ b ba 20 + h.c, 



(4) 



where A a = A^ cosh(r) — e l< ^Xi y sinh(r), A& = A 2j/ cosh(r) — e ll ^\ 2x sinh(r), and we have 
assumed Ai y cosh(r) — e l ^X\ x sinh(r) = A 2x cosh(r) — e l ^X 2y sinh(r) = 0, which implies that 



r = arctanh 
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(5) 



and 



<£>lz - (fly — - {^2x - ^2y) ■ 



(6) 



In this way, the squeezing factor r and the squeezing angle can be manipulated, respec- 
tively, by the intensities and the phases (fj a of the classical fields. When we take into account 
the atomic decay of levels |1) and |2), decay rates Ti and T 2 respectively, the dynamics of 
the system, in the transformed picture, is determined by the master equation 



P = - 



h 



E u p +C 1 p + C 2 p, 



(7) 
The 



where L x p = -f (2a 01 pa w - a u p - pan) and C 2 p = -f {^mP^ - V22P - PV22) 
steady state of Eq. (JTj) is the vacuum for both modes and |0) for the electronic state. 
We have assumed a strong decay of both excited electronic states once, as pointed in Ref. 
fis| . we need two distinct dissipation channels (£\p and £ 2 p) to protect a two-mode quan- 
tum state of a trapped ion. (Without this assumption we can not ensure that the steady 
state of both modes, in the transformed picture, is the vacuum state.) By applying the 
reverse unitary transformation, it is readily shown that the steady state of this system is 



p {t - 00) = S ab (OpSl(0 = S ab (0 |0, 0) (0, 0| 5^(0 ® |0> 



(8) 



which is a pure state for the vibrational modes a and b, i.e., exactly the two- mode squeezed 
vacuum state = S n tanh n (r) / cosh (r) |n, n) ab . The degree of squeezing r is determined 

and Ajq- = tTjj a Qj a . 
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by the amplitudes of the classical fields fl, a since tanh (r) 
This steady state does not depend on the initial electronic or motional state of the ion. 
Thus, the ion does not nave to be cooled to the fundamental state in order to prepare such 
a state. Also, as the entangled state is generated through the engineered reservoir, there is 
no requirement for a precisely timed interaction between the ion and the laser fields and the 
degree of entanglement (r) is determined only by the ratio of the amplitudes of the classical 
fields (see Eq. ((5])). In this scheme, the TMVS state is generated when the system reaches 
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261 ]. the time needed for the system to reach 



A, 



A, 



> T, this time will be 



the steady state. As pointed out in Refs. 
the steady state is defined by the atomic decay rate V. For 
of the order of a few times 1/r. In Ref. 27| Tang et al. showed how to generate the TMVS 
state in a single ion in a two-dimension trap inside a non-ideal optical cavity. Differently of 
our scheme, where the required dissipation channels are played by the decay of the excited 
electronic levels, in Ref. [271] the required dissipation channel is played by the decay of the 
cavity mode. 

To check our result we solve numerically the master equation for our system in the 
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interaction picture 



i 

P=~T i H I' P\ + + ^2P + £abP, (9) 



where Hi is given by Eq. (jHJ). In this equation we have introduced the Liouvillian L a bp 
which describes the action of the thermal reservoir on both atomic motions, i. e., 

ot=a,b ^ 

_l__t^Ta (2cJpa — aa^p — paa^) j , (10) 

n t /j being the mean number of quanta of the thermal reservoir and 7 a (7^) the decay rate 
of the vibrational mode x (y). We start with the ion in the internal ground state |0) and 
both modes in the thermal state, p a &(0) = p a <S> p&, with p a = Pb = S.^ , |n) (n|, n 

being the initial mean number of quanta for each mode. To solve numerically the master 
equation (j5J) we adopt the same coupling, X lx = X 2y = A, X ly = X 2x = Atanh(r), and the 
same decay rate for both excited electronic states, I\ = T 2 = T, and the same decay rate 
for the vibrational modes, 7 a = 7ft = 7- In Fig- 2, we have plotted the mean number of 
quanta of mode a against time (the evolution of mode b is identical) for different values of 7 
and for n = 2, n t h = 0.5, T = 10, and r = 1 (which implies tanh (r) = ^ = ^ ~ 0.76). 
For an ideal two-mode vacuum squeezed state, the mean number of quanta for each mode 
is (n a ,b) — i^tanh^ ' wn i cn for r = 1 gives (n a ^) ~ 1.4. We can see in Fig. 2 that, 
for 7 = 0.001A and 7 = 0.01A, a mean number of quanta close to this value is reached 
asymptotically, but this is not so for 7 = 0.1 A, because of the competition between the 
engineered and thermal reservoirs. Another parameter we have used to analyze the fidelity 
of the generated state is the total variance ({Au) 2 + {Av) 2 ) of a pair of EPR-like operators 
u = \e\ x a + \xb and v = |e| p„ — hPb |30] . with x a = (a + / \/2 and p a = —i (a — aty / \/2, 
a = a, b. According to Ref. [30], a two-mode Gaussian state is entangled if and only if 
((Am) 2 + (Av) 2 ) < e 2 + 1/e 2 . For e = 1 and an ideal two- mode vacuum squeezed state, the 
total variance is ((An) + (Av) ) = 2e , which for r = 1 gives us ((Am) + (Av) ) ~ 0.27. 
As we can see in Fig. 3, this value is reached and approximately reached for 7 = 0.001A 
and 7 = 0.01A, respectively. Again, for 7 = 0.1A, the action of the thermal reservoir does 
not allow the ideal generation of the two-mode entangled state. Instead of applying unitary 
transformations to the density matrix p, which led to Eq. (j3J), and thus making it easy to 
find the steady state, we could have proceeded by looking for an engineered Liouvillian for 



the engineered reservoir, as in Ref. [19j]. For an atomic decay rate T much stronger than 
the effective coupling A and the decay rate of the vibrational modes 7, the effective decay 
rate for the engineered reservoir is given by T eng = 4A 2 /r. In our numerical solution of the 
master equation ([9]) we have used T = 10A, which results in T eng = 0.4A, that is of the same 
order of magnitude as 7 = 0.1A. Then, for this value of 7, the influence on the generated 
TMVS state of the natural reservoir is almost the same as that of the engineered reservoir. 
Hence, to minimize the influence of the natural reservoirs, we must have T eng = 4A 2 /T ^> 7. 
For example, for 77 = 0.1 (Lamb-Dicke limit), Vt ia ~ 1 MHz, and T ~ 1 MHz, which can 
easily be achieved with current technology, we have |A| ~ 0.1 MHz and T eng = 40 KHz, 
which is much stronger than 7 ~ 2 KHz, found in present-day experiments. (The chosen 
values above also satisfy the requirements for the approximations employed to obtain the 
effective Hamiltonian: for u x ~ u y ~ 30 MHz >> \rjQ\ ~ 0.1 MHz and T ~ 1 MHz.) 

Summarizing, we have presented a simple scheme to prepare a two-mode vacuum squeezed 
state for the 2D motion of a trapped ion via the generation of an artificial reservoir. Our 
scheme is robust against the usual mechanism of dissipation and could be implemented with 
the present-day technology and we hope it could be employed to test experimentally the 
entanglement properties of Gaussian states. The approximations employed in this work were 
validated by numerical calculations, which showed complete agreement with the analytical 
solutions. To prove the engineered two-mode state a tomographic method couldbe employed 
that enables the Wigner function of the entangled state to be reconstructed |31| . 
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No. 2005/04105-5), and Brazilian Millennium Institute for Quantum Information. 
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Figure Captions: 

Fig. 1: Atomic levels of the trapped ion. The ground state |0) is coupled to the excited 
states |1) and |2) through laser fields. 

Fig. 2: The time evolution of the mean number of quanta, (a^a), of the vibrational mode 
x for T = 10A, n t h = 0.5, r = 1, and three values of the decay rate of the vibrational modes: 
7 = 0.001A (solid line), 7 = 0.01A (dotted line), and 7 = 0.1A (dashed-dotted line). The 
dashed line (straight line) represents the expected value. 

Fig. 3: The time evolution of the total variance ((Am) 2 + (Av) 2 ) for the same parameters 
used in Fig. 2. 
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